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Abstract 

We calculate the cross section of the process e'’'e“ —> ti near threshold by re¬ 
summing Coulomb-like terms with next-to-next-to-leading (NNLO) accuracy. The 
nonrelativistic Green function formalism and the method of “direct matching” are 
used. The NNLO correction turns out to be large, of the same size as the NLO 
correction. It changes the position and the normalization of the IS"—peak. The 
obtained results are compared with results existing in the literature. 
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1. Introduction. 


A detailed study of the process —>• ti will be performed at the Next Linear 

Collider []^. As a result of this study, precision measurements of the top quark mass m 
and width L and of the QCD coupling constant as{m) can be done at NLC. 

The cross section of the top quark production at LO and NLO was studied in detail 
in 0, B g. The general approach for the calculation of the cross section e'''e —> tt 
was developed in |^, which pointed out that the large width of the top quark serves as a 
cutoff for long-distance effects in this problem. This allows for reliable predictions for the 
cross section to be obtained by using perturbative QCD. However, it is well known that 
the standard QCD perturbation theory does not work in the threshold region. Even at 
LO approximation one has to resum all (ag/v)^ terms. The result of such resummation 
is proportional to the nonrelativistic Coulomb Green function at the origin. At NLO and 
NNLO one has to resum all (as/n)^(l, a^,, n, terms in order to calculate the 

cross section in the threshold region. 

The approach to the NNLO calculations was suggested in ^ and consists of two steps. 
First, one calculates the cross section in NNLO order far from the threshold, in the region 
Os -C 'y -C 1, ( n = — Am?/s is the relative velocity of top quark and antiquark), 

where the resummation of the Coulomb singularities is not needed. Then one considers 
the nonrelativistic cross section using an infrared cutoff and matches it with the cross 
section in the region Os -C n -C 1. 

The Abelian part of the corrections to the cross section near threshold proportional to 
the structures Cp and CpTNi were calculated in [^] and 0, respectively. The complete 
NNLO correction to the cross section for the threshold in the region as -C n <C 1 was 
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calculated in | 
threshold was considered recently in two papers [pA 


approaches. The difference of the numerical results between |T2| and is expected to 
be small, but actually the difference is large and amounts to 5 — 15%. 

The main purpose of this letter is to calculate the cross section of the process ti 

in the energy region close to the threshold by resumming Coulomb-like terms with next- 
to-next-to-leading accuracy. 

We use two alternative approaches and clarify the origin of the difference between [0 
Our hnal numerical result differs from both calculations ([1^] and |T^). 


and 

In section 2, we discuss basic elements needed to calculate the NNLO cross section. 
We consider the matching procedure and hx the short distance coefficient in section 3. In 
section 4, we present the hnal formula for the cross section and results of the numerical 
analysis. 


2. Basic elements of Rnnlo- 


The cross section of the process e’''e“ —ti in the near threshold region can be obtained 
by using the optical theorem and an expansion of the vector current The cross section 
reads 


247r r 

R = a(e+e" ^ ti)/(Tpt = G(ro)Im[(l - ^^)G(ro,ro|E -F iT) 
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where apt = 47 rQ!^/ 3 s, eg is the electric charge of the top quark, is the number of colors, 
y/s = 2m + E is the total energy of the quark-antiquark system, m is the top quark pole 
mass and F is the top quark width. G{f,r'\E + iF) is the nonrelativistic Green function, 
which satishes to the Schrodinger equation 

{H - E-iT)G{f,f'\E+ iT) = (2) 

The Green function is divergent at the origin r —0 because the potential contains 1/r^ 
terms. These singularities are regularized and factorized into the short distance coefficient 
C'(ro). 

In eq.(l) we have already used a complex energy, E + iF, in the Green function, ac¬ 
cording to 0 . This structure appears in the top quark propagator through the top quark 
self energy at LO, NLO and NNLO. The non-factorizable corrections were studied at 
NLO in |T^, 1^ and it was shown that they cancel in the total cross section. It was 
argued in |]T^ that the cancellation appears in all orders; therefore we do not consider 
nonfactorizable corrections here. 


2a. Hamiltonian and QCD potential. 


The nonrelativistic Hamiltonian of the heavy quark-antiquark system reads 

H = Ho + lT(r), Ho = - + V{r), (3) 

m 

here V (r) is a static QGD potential of the heavy quark-antiquark system at NNLO order: 

V{r) = ^^(2/?oln(/iir)-F ai) (4) 

2 

+ (^) (/^o(41n^(hi^) +-^) + 2(/3i-I-2/3oai) ln(/iir)-I- 02 )), 


here pi = f^s is the normalization scale, 7 is the Euler constant and ag = as{^s). 
The coefficients Oi and 02 were calculated in |^, |^, respectively 
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TrNl )' - (f - I6C3) CrTrNr. 


(5) 


The hrst two coefficients in the expansion of the QGD /d-function are 
11 4 34 20 

Po = -J^A - -NlTr, a = - -GaTrNr - AGrTrNr. ( 6 ) 

The color factors are Gr = 4/3, = 3 and Tr = 1/2. Nr = 5 is the number of light 

quarks. 
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2b. Breit-Fermi Hamiltonian. 


The function W (r) in eq.(^) is the QCD generalization of the QED Breit-Fermi Hamil¬ 
tonian [00- We consider here quark-antiquark production in the S'—wave mode. The 
Breit-Fermi Hamiltonian for the hnal state with L = 0, S*^ = 2 reads [^, 


W{r) = 


—H 


Am? 


llTrCi^Os Cpa^, 

H- ^^—o(r) - 




2m^ 
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CaCfciI 

2mr‘^ 


(7) 


Let us demonstrate how the Breit-Fermi Hamiltonian can be reduced to a convenient form 
for the numerical and analytical evaluation. First we note that eq. (@) can be rewritten 
in the following way 


W{r) 




llirCpas .5 CA.V^ir) 


( 8 ) 


Here we used Fo(r) = _ Hes® xhe hrst and second terms in eq.(||) can be simplihed by 
using equation of motion, [H — E)G{r', r"\E) = 6{r' — r"), E = E + ^T, but the third 
and fourth terms are related by the commutation relation 


[Ho, iPr] 


47r6(r) CpO-s 
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with ipr = Expressing the term with 1/r^ through 6(f) function and using the 

commutation relations, we have (assuming r ^ 0) 
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NNLO 


pr\E) = + 


( 10 ) 
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G^°(r,r|F) 
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We have omitted some unimportant surface terms above. Here Gc(T,r\E) is the Green 
function calculated only with the Coulomb potential and without the Breit-Fermi Hamil¬ 
tonian. The LO Green function can be written in the form iv = \ —) 


\ rri'' 


G^^( 0 , 0 \E) = ^(iv - GFas(\og(-^) -1-7^ + T(1 - (11) 

where /ij is a factorization scale, which disappears in R after performing a matching 
procedure. Employing this form of the LO Green function, we get the NNLO Green 
function in the form 
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NNLO 


(0,0|F) = 


riNNLO 

Lrc 


(0,0|B)-G7(0,0|B) 


( 12 ) 
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Equations ([l0|) and (^) are in agreement with |T^ but they are obtained in a different 
way. Expanding eq. m in vas, we have 


^iV7VLO(o^O|E) 

5(^)1 


r^NNLO 

Urc 


(0,0|E) + 5(G)i + 5(G)2 


(13) 
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These expressions are our analytical results for the correction to the Green function at 
the origin. 

An alternative strategy, more consistent one (adopted also in |^]) is to keep the Breit- 
Fermi correction together with the QCD potential in Schrodinger equation and to solve 
the Schrodinger equation numerically. The main reason for doing this is because the 
expansion procedure is not safe for energy denominators of the Green function 


G(r,r|E)=^ 


IV’n(r)P 


n 


E-iV 


+ 


dk \^|Jk{r)\‘^ 

2tt Ek — E — iT 


(14) 


especially if the energy is close to the energy of resonant state. 

Therefore, we re-express the 6{r) function in eq.(^ by the term with 1/r^ by using 
the commutation relation (^ 




,2 ty.TW 

^3 Cf 2m ■ 


(15) 


We consider the first correction to the Green function G{r',r"\E) originated from lT(r) 
Using the equation of motion we obtain 

' dfG{r',r\E)W{r)G{r, r"|E) = r"\E) + (16) 


drG{r\r\E)[ - 

4m 


2 m^ mr' ' '' ' 12m 

E^ 3VqE , ,2 , GA.V^ir) 


2m ^3 Gf^ 2m 


]G(r, r"\E). 


The surface term with ip^ does not coincide with the one from [|^. We see from the 
second line of the eq. (pTl) that the problem of evaluating the Breit-Fermi correction is 
reduced to solving the following equation 


{Hi-E^)Gi{r,r'\Ei) = 5%r-r'), 


with new energy Ei = E + ^ and with new Hamiltonian 


4m 


= - + U(r) 
m 


3VoE 2 GA.Vp^r) 
2m 3 Gf 2m 


(17) 


( 18 ) 


The numerical solution of (P!7[) gives us the Green function at NNLO order. 
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3. Matching procedure. 


Now we consider the matching procedure in order to £x the short distance coeffi¬ 
cient C{r) in eq.(|^. First we calculate the hrst-order correction to the Green function 
originated from W{r). This is given by the following integral 


- J d^rGoiro,r\E)W{r)Goir,ro\E) 

with the free Green function (with k = mv) 

m /sin{kr) , , sin{kr' 


Anrr 


Go{ry\E) = 0(r' - r) -f - r')). 


k 


k 


(19) 


( 20 ) 


Then we compare the result of this integration with the NNLO QGD radiative correction 
to the cross section e’''e“ —>■ tt in the region -C n -C 1 |]^. Iterations of the Goulomb 
potential give only terms proportional to 1/n or 1/n^ at NNLO in this kinematical region. 
These terms are trivially identihed in R. We obtain the following short distance correction 


G(r) = 1 - 


iGpOisifik) 


+ 


C2W( 


GF<^s{k'h)G 


). 


with 


TT ^ TT 

G2(r) = log(r/a) A 2 log(m//Uh) + A3 


2 Gjp p2-7B 

A,=y(GA + ^), « = A2 = 2/3o, 


and with 


A3 = GfG^^ + GaG^ + Tn^G^^ + TnE^G. 
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Here is the hard normalization scale. We should note that matching of the analytical 


expressions of the eqs. (pTD - (pAD gives similar results but with substitution r/a —>• m/iif. 


These coefficients are in agreement with results obtained in [1^, 

4. Numerical results. 

The hnal expression for the NNLO cross section is given by the following equation 


Rnnlo{E) = (i-^^^ + G, 

^ TT 


Cf<^s\‘^\ Stt 


- ^)Gi{ro,ro\E^)), (26) 
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with G'i(ro,rol-Ei) being solution of the eqs. ([T7|) - (p!8D and with C2(ro) from 
The analytical result for Rnnlo is given by the same formula (P^D with G(0,0|^) from 
eqs.(p!^ ) -([T^ ) and with substitutions: r —> ma/fXf in C{r) and (1 — |^) —>• (1 — |^) in 

eq.®. 

In order to solve the Shrodinger equation numerically we have written two programs in 
FORTRAN and MATHEMATICA by using numerical schemes described in Q and |T^ . 


We have checked that the programs reproduce the results of the analytical expression for 
the Green function in the pure Coulomb problem and numerical values for NLO cross 
section from B3I3. The results of both programs for the NNLO correction give 


the same result; therefore we are confident in our numerical results. 

The Fig.l shows our final results for Rmnlo{E) as a function of the nonrelativistic 
energy E = y/s — 2m. We compare the NNLO results with LO and NLO curves, for the 
soft scale = 50, 75,100 GeV. We have chosen m = 175GeV, F = 1.43GeR, tq = a and 
as{Mz) = 0.118. We see that NNLO correction is of the order 20% and as large as NLO 
one. The IS* peak is shifted towards smaller energies. The dependence of the NNLO cross 
section on the parameter is weaker than in LO, but more robust than in NLO. 

The dependences on the factorization scale tq and hard normalization scale are 
much smaller than the dependence on /r^. We demonstrate that fact in Fig 2., where we 
plot Rnnlo{E) as a function of energy at rg = 2a, a, a/2. In Fig.3 we show Rnnlo{E) at 
fih = ‘2m, m, m/2. In Fig. 4 we demonstrate Rnnlo{E) at different as{Mz) = 0.116, 0.118, 0.12. 
We see that Rnnlo{E) is very sensitive to the value of the QCD coupling as{Mz). Gom- 
paring our results (Fig.l) with the numerical results of [|T^] we have found that they differ 
by about 4 — 6% at the IS* peak and above it and by 10 — 20% bellow the IS' peak, where 
the cross section is quit small and NNLO correction is large. The difference of our results 
and the results of 


131 appears because of error in numerical solution of the Schrodinger 


equation. 

If we expand the Green function in the Breit-Fermi Hamiltonian we obtain the result 
that is in agreement with result presented in |^, see eqs. (|T^-(0). However, the 


expansion of the energy denominator of the Green function is not correct in the resonance 
energy region, where an expansion parameter, pf'" , is large. We do not use it in our 


calculation, and prefer to solve Schrodinger equation numerically with effective potential 
contained the Goulomb potential and the Breit-Fermi potential. Therefore, our hnal result 
differs numerically from both results presented^ in [|I^ and . The qualitative conclusion 
about the large size and the positive sign of the NNLO correction is in agreement with 
and ||I^ . 


5. Conclusion and outlook. 

In conclusion, we have calculated the total cross section, resumming Goulomb-Iike terms 
with next-to-next-to-leading accuracy. We used the nonrelativistic Green function formal¬ 
ism and the method of “direct matching”. The NNLO correction turns out to be large, of 
the same order as the NLO correction. It shifts the position of the IS*—peak and changes 
the its normalization. 


The comparison of obtained results with the results published in [12, 13| shows that 


^ After this paper has been accepted for publication, the authors of have agreed with the numerical 
results obtained in this paper. 
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our final results are in qualitative agreement, but differ numerically from that of 

0 . 


and 


Finally, let us mention some open questions in this held. It is quite important to 
consider NNLO correction to the total cross section in the scheme with the running QCD 
coupling as and in the scheme with the low-scale running top quark mass, for example 
mpsin) 0. It would be interesting in future to examine the size of the NNLO corrections 


to the differential distributions, for example to the distribution over spatial momentum. 
It is necessary to analyze nonfactorizable corrections to the differential cross section 
at NNLO. In view of the fact that NNLO correction is large it would be useful to check 
NNLO correction to the static QCD potential, the coefficient 02 
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Figure 1: i?(e+e“ —*■ ti) for the LO (dashed-dotted lines ), NLO (dashed lines), NNLO 
(solid lines) approximation as a function of energy E = y/s — 2m, GeV. In all cases we 
use mt =175 GeVjF^ =1.43 GeV, as{mz) =0.118 but different values of the soft scale 
Us = 50, 75,100 GeV . 
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Figure 2: R{e^e~ tt) at NNLO at different values of factorization scale tq = 2a(dashed- 
dotted line), tq = a (dashed line) and ro = a/2 (solid line). We use rrit =175 GeV, F* =1.43 
GeV, Us = 75 GeV,as{mz) =0.118. 
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Figure 3: -R(e’''e“ — ti) at NNLO at different value of the hard normalization scale 
Hh = 2m(dashed-dotted line), /ih = m (dashed line) and Hh = raj2 (solid line). We use 
mt =175 GeV, F* =1.43 GeV, /is = 75 GeV,as(m^) =0.118. 
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Figure 4: R{e^e~ —> ti) at NNLO at different values of the QCD coupling constant 
(y-sijnz) =0.116 (solid line), asijnz) =0.118 (dashed line),as(m^) =0.120 (dashed-dotted 
line). We use mt =175 GeV, F* =1.43 GeV and = 75 GeV. 
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